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Abstract. We present here a new determination of the ratio of energies radiated by active galactic nuclei and by 
stars and discuss the reasons for the apparently conflicting results found in previous studies. We conclude that 
the energy radiated by accretion processes onto super massive black holes is about 1 to 5% of the energy radiated 
by stars. We also estimate that the total mass accreted on average by a super massive black hole at the centre of 
a typical 10 11 Mq galaxy is of about 710 7 M©. 
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1. Introduction 

The universe is composed of various objects with a wide 
variety of emission processes. However, photons have three 
main origins: the big bang as the source of the cosmic 
microwave background (CMB), nuclear reactions in stars 
and gravity through accretion onto compact objects, es- 
pecially super massive black holes (SMBH) . Currently (at 
z ~ 0), the ratio of the number of photons of CMB and 
stellar origin in the universe is of the order of 400. It is 
more difficult to estimate the ratio of the nuclear energy 
radiated by stars and the gravitational energy radiated 
by active galactic nuclei (AGN). Several recent studies 
have shown that the energy r eleased by nucleosynthesis 
is evolving with cosmic time llMadau et all Il996l [l998: 
iTrentham et alJll999t ISomerville et alJl200l|) . It further 
seems that the star formation history roughly matches the 
evolut ion of energy release through accretion onto black 
holes !|DunlorJll998t iFranceschini e t al. 1999), since both 
processes peak around z ~ 2. It is therefore possible to 
relate both phenomena and to estimate the energy ratio 
radiated by AGN and by stars over the history of the uni- 
verse. 

The determination of thi s ratio has bee n addressed 
several times in the literature. iDunlod l)l998|) constructed 
two models of the radio luminosity function of a sample 
of radio-loud quasars (RLQ), the first only considering lu- 
minosity evolution and the second combining luminosity 
and density evolution. He then related this luminosity to 
the mass accreted onto the central black hole and corn- 
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pared these curves with the star formation history. The 
two curves appear to be correlated, suggesting that when 1 
M Q is accreted by a SMBH in a radio loud quasar, 10 7 M© 
are use d in the sta r forma tion process. Based on this re- 
lation, ICourvoisierl l|200lj) derived a ratio of the energy 
radiated by RLQ and stars of -Erlq/ 'E Bt ^ B ~ 10 -5 by as- 
suming an accretion efficiency of 10 % for the RLQ and a 
stellar population made of 10 solar mass stars each radiat- 
ing 4 1 52 ergs over the ir lifetime. Rather t h an co nsidering 
only radio loud AGN, IFranceschini et all 1^399) use the 
0.5 — 2 keV X-ray emission as a measure of the energy ra- 
diated by AGN. They find that when 1 M Q yr -1 Mpc" 3 is 
absorbed for star formation, the 0.5 — 2 keV volume emis- 
sivity from AGN is 2.4 1 40 ergs s _1 Mpc -3 . Considering 
that type I and II AGN bolometric luminosity is 250 times 
the 0.5 — 2 keV luminosity and a stellar radiative effi- 
ciency of 0.001, IFranceschini et alJ l|1999|) finally obtain 

i bol (AGN)/i bol (SFR) ~ 0.1. 

A different appr oach was followed by 
iFabian k. Iwasawal l)l999|) . Their estimate of this ra- 
tio is based on relations linking the bulge mass of a 
galaxy to both the mass of its central black hole and of 
its stars. They obtained a value of ~ 0.2 for the ratio of 
the energy radiated by accretion processes and by stars, 
when assuming an accretion efficiency of 10 % for the 
AGN and the fact that one tenth of the stellar mass is 
used for nu clear fusion with a n efficiency of 0.6 %. 

Finally, lElvis et a l. ( 2002) estimated that AGN con- 
tribute by at least 7 % to the total luminosity of the uni- 
verse as derived from the diffuse background from sub- 
millimeter to ultraviolet wavelengths. They obtained this 
result by first deriving a lower limit of the AGN X-ray 
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emission from the X-ray background (XRB) and by apply- 
ing a bolometric correction determined with the average 
spectral energy distribution of quasars. 

The aim of this study is to compare the various studies 
described above and to derive a new estimate of the ra- 
tio of gravitational energy released around SMBH to the 
energy released by nuclear fusion in stars. We first derive 
the radiation energy density due to stars in Sect. El then 
the corresponding value for accretion by SMBH based on 
the XRB in Sect. |31 The obtained ratio is compared with 
previous studies in Sect. where we try to identify the 
origin of the discrepancies. 

2. Energy radiated by stars 

To estimate the energy radiated by stars in the universe 
we need to know both the energy release of a typical stellar 
population and the overall star formation history. Below, 
we start with the calculation of the stellar energy release, 
while the evolution of the star density will be described in 
Sect. IP 

2.1. Stellar energy release 

We use the Starburst 99 l models of lLeitherer et all 1^999) 
to determine the typical stellar energy release. These mod- 
els predict the spectrophotometric evolution of starburst 
galaxies between 10 6 and 10 9 years after the onset of 
star formation bas ed on the stellar evolution models of 
the Geneva group llSchaller et al.lll992HCharbonnel et alJ 
Il993t ISchaerer et alJll993albt iMevnet et all 11994^ . They 
consider the atmosphere mo d els of iLeieune et aP l|l997l) 
and those of ISchmutz et all l|l992h when the mass loss 
becomes important. A simple black body is used for 
cool stars with additional nebular continuum including 
free-free interactions below 912 A and free-bound interac- 
tions above. These models have been computed with the 
isochrone synthesis method and are optimized for massive 
stars. 

Since we are only interested in the total energy release 
of a typical population of stars during its whole life, we 
only consider the inst antaneous star formation models of 
iLeitherer et al.l £!)!)!>! because in this case most of the 
energy is released before 10 9 years. To assess the effect of 
changing the powerlaw index a of the initial mass func- 
tion (IMF) of the stars we consider both a Salpeter IMF 
(a = 2.35) and a steeper Scalo IMF (a — 3.3). In both 
cases, the cutoff masses are chosen as M low = 1 M Q and 
M up = 100M Q . The effect of changing the metallicity Z is 
taken into account by considering four different metallic- 
ities: Z = 0.040, 0.020(= Z & ), 0.008 and 0.001, but with- 
out chemical evolution in the models. As an example, we 
show in Fig. ^ the evolution of the bolometric luminos- 
ity of a 10 6 M Q star cluster formed instantaneously with 
a solar metallic i ty [Z = 0.020) according to the model of 
ILeitherer et al.l l)l999j) . We extrapolate the bolometric lu- 

1 http:/ /www. stsci.edu/science/starburst99/ 
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Fig. 1. Evolution of the bolometric luminosity of a 
10 6 M Q star cluster formed instantaneousl y with a metal- 
licity Z = 0.020 according to the model of lLeitherer et alJ 
l|l999l) . The solid line is for a Salpeter IMF (a = 2.35) and 
the dashed line is for a Scalo IMF (a — 3.3). Only stars 
in the mass range between 1 and 100 M are considered. 
Beyond 10 9 years (dotted line), the model data have been 
linearly extrapolated up to 10 10 years. 

Table 1. Energy radiated per unit of solar mass by stars 
(£ rad /M ) atars and the corresponding energy density (7 stars 
for different metallicities Z and for a Salpeter (a = 2.35) 
or a Scalo (a — 3.3) IMF. The last column shows the cor- 
responding energy density ratio radiated by AGN and by 
stars. 



stars £4tars UAG~N/U B t ars 

10 51 crgsM Q _1 10 60 crgsMpc" 3 (10~ 2 ) 

0.001 235 4~16 L29 2l 

3.3 4.69 3.33 1.1 

0.008 2.35 3.46 1.07 3.4 

3.3 3.39 2.40 1.5 

0.020 2.35 3.09 0.96 3.8 

3.3 2.92 2.07 1.8 

0.040 2.35 2.33 0.72 5.1 

3.3 1.71 1.21 3.0 



minosity from 10 9 to 10 10 years with a power-law in order 
to include the energy radiated during the final stages of 
stellar activity. This ext rapolation is in good agr eement 
with the earlier study of lCharlot fc Bruzuall (|l99l|). 

The total energy radiated by a 10 6 M Q star cluster 
is obtained by integrating the bolometric luminosity over 
time. The results are presented in Table \l\ior the instan- 
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taneous stellar formation law, four different metallicities, 
and two different IMF. Since the bolometric luminosity 
is dominated by the energy radiated by the most mas- 
sive stars, the flatter powerlaw index a of the Salpeter 
IMF provides a higher energy release than the Scalo IMF 
until oj 10 s years. Afterwards, the radiated energy being 
provided by less massive stars, a Scalo IMF gives more 
energy. Globally, the difference between the total energy 
radiated by a cluster of stars for a Scalo or a Salpeter IMF 
is quite small. On the other hand, the effect of increasing 
the metallicity is to decrease the stellar energy release. 
This effect can be understood as being due to an incre ase 
of stellar opacity with metallicity llMowlayi et al Il998h . A 
quick comparison with ICourvoisierl l)200l|) shows that the 
values in Table n are consistent with his rough estimation 
that stars radiate ~ 10 51 ergsMg 1 over their entire lives. 



2.2. Evolution of the stellar density 

iMadau et al.l l|l996|) first computed the star formation as 
a function of redshift. Since then, several authors have 
added new points to his diagram of the star formation 
rate (SFR) per unit of comoving volume as a function of 
the redshift. Different tracers have been used to derive the 
SFR, but now, it is common to use galaxy luminosities at 
different wavelengths. They can be converted into stellar 
formation rates using stellar population and galaxy spec- 
tral models, stellar formation scenarios and various IMF. 
IMadau et all l|l998(l propose conversion factors from lu- 
minosity to star formation rates (L vv = C ■ SFR) with 
differ ent values of C f or a Salpeter IMF and a Scalo 
IMF. ISomerville et al.1 (|200ljl compiled all the observa- 
tions made in this way and present a homogeneous table 
of the comoving SFR density data for different cosmologi- 
cal models and a Salpeter IMF. Fig.[2]shows those data in 
the case of an Einstein-de Sitter cosmology and a Salpeter 
IMF. We calculate the stellar density p stars by integrating 
the star formation rate /9 stars over the whole cosmic time: 



,dt. 



(1) 



If we consider a Friedman cosmological model, a fiat uni- 
verse, q = 0.5, and H = 50 kms _1 Mpc -1 , Eq. QJ be- 
comes: 



Psta 



Astars (Z) 



dz 



H (l + z) 



2.5 



(2) 



With z max = 5.5 and p s tars(z) being the solid line adjusted 
to the data points of Fig. we obtain: 

p stars = 3.1-I0 8 M Q Mp C - 3 . (3) 

This value is an average over the st ar formation history . 
Using the conve rsion factors given bvlMadau et alJ l)l998(l 
and the table of lSomerville et al.l(|200l[) . we deduce a dif- 
ferent stellar density if we consider a Scalo IMF. Thus, 
for a Scalo IMF, an Einstein-de Sitter cosmology and 
Hq = 50 kms -1 Mpc -1 , we get: 

p stars = 7.1-I0 8 M o Mpc~ 3 (4) 



Z, -2 
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• Madau et al (1998) 
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Fig. 2. Evolution of the comoving SFR density for an 
Einstein-de Sitter c osmology and a Salpete r IMF from the 
data in Table A2 of lSomerville et alJ (|200lh . The solid line 
is a cubic spline matching as well as possible the observa- 
tional constraints. 

which is about twice the value for the Salpeter IMF. These 
results are only very weakly dependent on the cosmolog- 
ical model. Indeed, if we consider a ACDM cosmological 
model (fio = 0.3, il\ = 0.7) and calculate p s t ars , we find a 
difference that is negligible (~ 0.2 %) for both a Salpeter 
and a Scalo IMF compared to our previous values. 

2.3. Energy density of the stellar radiation field 

Having obtained the energy radiated by stars (see Tabled 
column 2) and the star density for a Salpeter IMF (see 
Eq. 01 and for a Scalo IMF (see Eq. 0J, we calculate the 
energy density due to stars as 

-E ra d 



^stars Psta 



Mr. 



ergs Mpc 



(5) 



The obtained values are reported in the fourth column of 
Tabled Since U stais is proportional to p sta _ IS and E la _ d (stars) , 
we note that U ata . ra decreases when Z increases following 
the E rai (stars) behavior and increases when the IMF slope 
is steeper. 

3. Energy radiated by active galactic nuclei 

To estimate the energy radiated by AGN, we use the 
observed X-ray background (XRB). This is motivated 
by the growing evidence that the XRB is emitted by 
discrete sources that are mainly AGN as proposed by 
ISetti k, Woltierl ljl989j) . By adding the contribution of 
type I and II AGN, the observed XRB spectrum can in- 
deed be well reproduced l|Comastri et alJll99.4 iGilli et ail 
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Il999i l200l|) . Furthermore, it seems that the contribution 
of both massive X-ray binaries and superno vae to the 
XRB is neglig ible l|Nataraian fe AhnainiK OOO). The emis- 
sion from the hot interstellar m edium is also small c om- 
pared to the AGN contribution l|Comastri et al.lll995l) . It 
is therefore possible to calculate the energy density in the 
AGN radiation field, CTaon, from the observed XRB spec- 
trum. iFab ian & Iwasawa (1999) give the following ana- 
lytic parametrisation that describes well the spectral en- 
ergy distribution of the XRB as observed by HEAO-2 and 
ASCA for soft and hard X-rays: 

1(E) = 9 • E-° A • exp(-£:/50kev) keVcuT 2 s" 1 sr" 1 keV _1 (^! 

To relate this spectrum to the bolometric intensity of 
AGN, we use the mean ratio from the 2 — 10 keV intensity 
to the bolometric intensity derived bv lFabian fe Iwasawal 

ET 



^[2-iokov]/^boi (RQQ) — 3.3%. 



(7) 



This ratio is based on a compilation of spec tral energy 
distrib utions of radio-quiet quasars (RQQ) bv lElvis et alJ 
l|l994h . which can be considered to be representative 
for the aver age emission of unobscur ed AGN. We fur- 
ther assume [Fabian fc Iwasawal l)l999j) 's hypothesis that 
the unobscured emission E 1(E) of AGN down to 2 
keV can be considered as constant with a value of 
38 keVcm~ 2 s _1 sr _1 which is the XRB intensity at 30 keV. 
We then calculate the 2 — 10 keV intensity of AGN as being 
61 keVcm~ 2 s _1 sr _1 . From this last value and Eq. J7J), we 
derive the bolometric intensity of AGN and we calculate 
the energy density of the AGN radiation field as 



problem; sometimes observational and sometimes approx- 
imative with the mass to energy conversion derived ac- 
cording to an efficiency e rad through E ra _ d = e rad Mc 2 . 

iFabian fc Iwasawal 1^)99) consider a typical galaxy 
and the contribution to the flux radiated both by the stars 
and the black hole hosted in this galaxy. They implicitly 
assume all galaxies with a black hole to host an AGN. 
However, the presence of a black hole does not mean that 
this galaxy is active. Some galaxies might have had an ac- 
tive phase, but are currently quiescent. Taking their ratio 
is equivalent to considering that every galaxy is currently 
tive and contributes to the radiated energy of AGN. 
is leads to overestimate the energy radiated by AGN 
and thus to an overestimation of the AGN-to-star radia- 
tion ratio. 

iFranceschini et afl ljl999j) tackles the problem from the 
observational point of view starting with the energy radi- 
ated by the stars and the AGN to find a link between 
both phenomena. This approach is motivated by the ap- 
parent similarity between the cosmic evolutio n of the SFR 
and the AGN activity. In ICourvoisierl l)200l[) . the sample 
of these active galaxies is based only on a radio survey. 
Therefore, this selection results in a sample of radio-loud 
galaxies which are only a small fraction of the whole AGN 
population. This leads to an underestimation of the accre- 
tion rate per unit of volume because there are many radio- 
quiet quasars that contribute to the energy radiated by ac- 
cretion processes that are not taken int o account. The use 
of the 0.5 — 2 keV volume emissivity bv IFranceschini et alJ 
199S 



AGN 



3.7 1 58 ergsMpc 



-3 



(8) 



radiate in the X-ravs jCcanastrietalJ 


Il995l 


[Gilh et alJ 


1999; 


Mivaii et alJ 


2000; 


iPompilio et al. 


2000; 


iGilli et alJ 


200l| while only including a negligible contribution of 



4. Discussion 



other objec 


like star clusters or massive X-ray binaries 


llGilli et alJ 


1999; 


iNataraian & Almaini 


2000; 


iGilli et alJ 



We have independently determined the energy densities 
emitted both by stars and AGN in the universe. The ob- 
tained ratio UAGN/U ataia is given in the last column of 
Table ^ for two IMF and four different metallicities. We 
note that all values are between 10 -2 and 510 -2 . They 
slightly depend on the IMF slope and on the metallic- 
ity. Moreover, the choice of the cosmological model has a 
negligible effect as mentioned at the end of section 12.21 
This ratio is expected to remain rather constant with cos- 
mic time because the star formation history is similar to 
the AGN luminosity evolution toward s higher redshifts 
l|Dunlodll998l: IFranceschini et al.lll999(l . 

4.1. Comparison between the different approaches 

Our results differ by several orders of magnitude from 
other estimations of the energy ratio radiated by AGN 
and by stars. The different r esults extend from 10~ 5 
llCourvoisierl l200l|) to IP" 1 llFabian fc Iwasawal Il999t 
IFranceschini et alJll999tlElvis et alJl2002|) through 4 10" 2 
(this work). We cannot directly compare these values, be- 
cause each study is based on a different approach to the 



120011). 

lElvis et"afl |2002) also consider the AGN emission 
based on the XRB. However, they compare it to the to- 
tal luminosity of the universe that is estimated from the 
diffuse background from submillimeter to UV wavelengths 
rather than only the stellar luminosity. 

We conclude that the various studies mentioned above 
actually measure different quantities. Therefore, it is not 
possible to compare directly the values obtained by differ- 
ent groups. 

4.2. Comparison between the different parameter 
values 

In addition to the differences in the approaches and the 
measured quantities pointed out above, various studies use 
different values for t he same parameter s . The stellar effi- 
ciency e rad used by 
but it is of 0.001 in 
the ratio of the 2 — 10 keV luminosi ty to the bolomet- 
ric lum inosity of the AGN is used in | F abian fc Iwasawal 
(1999) and in this work, while IFranceschini et al.l 1)19991) 



Fabian fc Iwasawal dl999l) is 0.0006, 

Similarly, 



Franceschini ct al. ( 1999 
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consider instead the 0.5 — 2 keV flux to derive a bolo- 
metric luminosity. Using the spectrum of Eq. © , we can 
convert their bolometric correction to the one based on the 
2-10 keV flux. Thus, we derive that th eir bolometric cor- 
rection differs from the 3.3% value of IFabian fc Iwasawal 
l)l999(l by a factor of 3. It seems therefore that there 
is an accumulation of different factors explaining the 
diverging results in the literature. By taking the val- 
ues of £ rar i = 0.0006 and Ln— iokeV]/-^boi = 3.3% from 
IFabian fc Iwasa wal <ll9 99 [) and repeating the calculation 



of IFranceschini et al J (fl99^ . we obtain a value of ~ 0.053 
for the ratio L bo i(AGN)/L bol (SFR), which tends to the value 
of 0.04 we obtain in this work for a solar metallicity and 
a Salpeter IMF. 

In order to expl icitly calcu l ate U. t ^ as defined here 
from the study of lElvis et alJ ^2002^ . we first subtract 
the quasar contribution to the total luminosity of the 
universe they give to get only the stellar background. 
We then obtain that 6.6 1 59 ergs Mpc -3 are radiated by 
stars which is quite similar to our result (see Table . 
Therefore, keeping their value of C/ A gn, we estimate the 
ratio U A GTst/U Bt!iIB ~ 0.067 that becomes ~ 0.046 if we use 
the energy radiated by stars obtained in our study for a 
Salpeter IMF and a solar metallicity (see Table instead 
of theirs. If we compare this to our results, we note that 
a difference also r esides in the value of C/agn ■ Considering 
lElvis et all l|2002|) 's values of 48 keVcm" 2 s" 1 sr" 1 at 30 
keV for the XRB instead of 38 keV cm~ 2 s _1 sr _1 and their 
correction of a factor 1.6 to the bolometric correction, we 
recalculate our value of C/agn applying the same method 
as seen previously in Sect. El We find t hat the AGN emit- 
ted energy density of lElvis et alJ (|2002j) is 1.5 times higher 
than ours. With this new derivation of C/ A gn, we calculate 
a ratio C/AGN/C/ ata r S ~ 0.031 for a Salpeter IMF and a so- 
lar metallicity that is a factor 2 lower than their value 
of 0.067. The difference comes from a lower C/ A gn and a 
higher C/ Stars , both effects combining to give a ratio twice 
lower. 

Finally, i t is w orth noting that recent observations by 
ISarzi et alJ (120011) lead to a smaller ratio of the central 
black h ole mass to the host bul ge mass than the value 
used bv lFabian fc Iwasawal l|l999l) . By using their new re- 
sul t of M BH /M h „i gf . ~ 0.002 the ratio -Eagn / E ettaa found 
by IFabian fc Iwasawal l)l999|) would have been of ~ 0.07 
instead of ~0.2. 



4.3. Relation between RLQ radio luminosity and AGN 
bolometric luminosity 



IFranceschini et ail (^399) and ICourvoislerl (j200lT) com- 
pared the star formation rate to the accretion history 
based on the observation either in the X-ray band or in 
the radio band. Every AGN is responsible for the X-ray 
emission but only a subset of AGN has a significant radio 
emission. Therefore, we easily get a link between these two 
classes of objects. When 1 M Q yr" 1 Mpc~ 3 is used by the 
star formation, 2.4 1 40 ergss -1 Mpc~ 3 are emitted by the 



us- 



AGN in the 0.5 - 2 keV band and 10 -7 M Q yr" 1 Mpc -3 
is accreted onto the central black hole of a RLQ. If we 
transform the 0.5 — 2 keV luminosity into the bolomet- 
ric luminosity using L[ .5-2kcV]/^boi(AGN) ~ 1.25% de- 
rived from Eqs ijfjj and Q and if we consider an accre- 
tion efficiency onto the black hole of 10 %, we obtain that 
3.4 10~ 4 Meyr^Mpc -3 is accreted by the AGN for a 
SFR of 1 M Q yr _1 Mpc~ 3 . Therefore, we can compare the 
luminosity of both the RLQ and the AGN as we get ~ 10 4 
of difference between them. When 1 ergs -1 Mpc~ 3 is ra- 
diated by the RLQ at 2.7 GHz, the complete population 
of AGN radiates about 3400 times more. 

Furthermore, w e can reconsider the method used by 
ICourvoisierl l|200lj) to determine U AGN /U ataTS . We have 
seen that we cannot directly compare his result to the 
others since they are not considering the same f amily of 
object s. Instead of using the analysis of RLQ bv iDunlonl 
(|l998fh we recalculate the rati o- of ICourvoisierl ll200 ll 
ing the analysis of the XRB of IFranceschini et al.l fa 
Therefore, we obtain a value of E AGJq / E BtaTB ~ 0.015 be- 
tween the energies radiated by AGN and stars instead of 
RLQ and stars. This value becomes ~ 0.020 if we use our 
result of the energy radiated by stars for a Salpeter IMF 
and a solar metallicity (see Table The initial inconsis- 
tency of a factor of a thousand has been reduced to only 
a factor of 2 between this last estimate and the one in 
Table □ 



4.4. Mass accreted by supermassive black holes 

Based on our previous results we can estimate the total 
mass accreted by SMBH. If we consider a galaxy with a 
mass of 10 11 M Q in stars, we can estimate the energy radi- 
ated by those stars from the value in Table^for a Salpeter 
IMF and a solar metallicity. Using the corresponding ra- 
tio of the AGN-to-stars energy release, we can then derive 
the energy radiated by the SMBH. By further assuming 
an accretion efficiency of 10%, we obtain that the total 
mass accreted by the SMBH is 6.5 10 7 M . In general, 
the typical mass accreted by a SMBH in a galaxy can be 
calculated according to 



M?g N (M Q ) 



- , 0.1 
2.2 • 10 7 I — 

M Gal 



a 



0.04 

(£ rad /M Q ) 



10 n M O/ / VlO 51 ergsM c . 



(9) 



We can also estimate this mass starting from the link 
given in the previous subsection. We have seen that when 
1 M Q of material is converted into stars, 3.4 10 -4 M Q is 
accreted onto SMBH. If we consider again a galaxy with 
a mass of 10 11 M Q in stars, we estimate that 3.4 10 7 M Q 
is accreted onto the SMBH that is half of our previous 
estimate. This difference by a factor of 2 is the same as 
for the ratio E AGN /E stllIS mentioned in Sect. 14.31 



6 



Zurita et al.: Ratio of energies radiated in the universe through accretive processes and nucleosynthesis 



5. Conclusion 

We derived the relative contribution of AGN and stars 
to the radiation energy density of the universe. The re- 
sults are given for different IMF and metallicities of the 
stellar population. The obtained values for the energy ra- 
tio released by AGN and stars are all between 0.01 and 
0.05. In the case of a Salpeter IMF and solar metallicity 
we obtain a ratio of 0.04. This result cannot be compared 
directly with previous studies because the approaches as 
well as the values of some parameters used in the calcu- 
lation differ from one study to the other. Actually, when 
using similar parameters and appropriate correction fac- 
tors it seems that all previous studies do converge to values 
between 0.02 and 0.07 for the ratio of energy radiated by 
AGN and stars. 

Since the CMB energy density is about ~ 0.2 eVcm~ 3 , 
the energy density due to stars is of the order of 
10 60 ergsMpc -3 - 0.02 cVcnT 3 and the AGN related 
radiation energy density is of about 8 10~ 4 cVcm~ 3 , the 
general picture resulting from this work is that the CMB 
contributes 10 times more than stars and 250 times more 
than AGN to the local energy density of background pho- 
tons. 

We also estimate that RLQ contribute about 3400 
times less than the whole AGN family to the total accre- 
tion power in the universe and that the cumulated mass 
accreted on average by a SMBH is of about 6.5 1O 7 M . 
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